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Possible physical realizations of the Tolman VII solution 
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The Tolman VII solution for a static perfect fluid sphere to the Einstein equations is reexamined, 
and a closed form class of equations of state (EOSs) is deduced for the first time. These EOSs 
allow further analysis to be carried out, leading to a viable model for compact stars with arbitrary 
boundary mass density to be obtained. Explicit application of causality conditions places further 
constraints on the model, and recent observations of masses and radii of neutron stars prove to 
be within the predictions of the model. The adiabatic index predicted is 7 > 2, but self-bound 
crust solutions are not excluded if we allow for higher polytropic indices in the crustal regions of 
the star. The solution is also shown to obey known stability criteria often used in modeling such 
stars. It is argued that this solution provides realistic limits on models of compact stars, maybe even 
independently of the type of EOS, since most of the EOSs usually considered do show a quadratic 
density falloff to first order, and this solution is the unique exact solution that has this property. 

PACS numbers: 04.20.Jb, 04.40.Dg, 95.30.Sf, 97.60. Jd 


I. INTRODUCTION 


The construction of exact analytic solutions to the Ein¬ 
stein equations has had a long history, nearly one hun¬ 
dred years to be more precise. However, in spite of the 
fact that the total number of solutions is large [1J and 
growing, only a small subset of those solutions can be 
thought of as having any physical relevance. Most solu¬ 
tions exhibit mathematical pathologies or violate simple 
principles of physics (energy conditions, causality, etc.) 
and are therefore not viable descriptions of any observ¬ 
able or potentially observable phenomena. 

Indeed, works that review exact solutions and their 
properties demonstrate the difficulties associated with 
constructing solutions that might be relevant to gravi¬ 
tating systems that actually exist in our Universe. Even 
in the simplest case of exact analytic solutions for static, 
spherically symmetric fluid spheres, it has been shown 
that less than ten percent of the many known solutions 
can be considered as describing a realistic, observable 
object. For example, |Delgaty and Lake using com¬ 
puter algebra methods, reviewed over 130 solutions and 
found that only nine could be classified as physically rele¬ 
vant [2]. A similar study by Finch and Skea arrived at the 
same conclusion [3|. The latter review also introduced an 
additional criterion that further reduced the number of 
physically relevant solutions to those that have exact an¬ 
alytic equations of state (EOSs) of the form p = p(p), 
where p is the fluid pressure and p is the matter density. 
This class of solutions was called “the set of interesting 
solutions.” 


In 1939 Tolman introduced a technique for construct¬ 
ing solutions to the static, spherically symmetric Einstein 
equations with material fluid sources [1]. That method 
led to eight exact analytic expressions for the metric func- 
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tions, the matter density and in some cases the fluid pres¬ 
sure. Beginning with an exact analytic solution for one of 
the two metric functions, an expression for the mass den¬ 
sity could be obtained by integration. With expressions 
for the density and the first metric function in hand, an 
analytic expression for the second metric function could 
be obtained. This often required an appropriate change 
of the radial variable to obtain a simple integral. All 
functions could then be written as explicit functions of 
the radial coordinate r. While the fluid pressure could, in 
principle, be obtained from the metric and density func¬ 
tions, Tolman chose not to evaluate the fluid pressure 
in some cases due to the fact that to do so would lead 
to mathematically rather complicated expressions that 
might be difficult to interpret. 

Of the eight solutions presented in his paper, 
three were already known (the Einstein universe, the 
Schwarzschild-de Sitter solution, and the Schwarzschilcl 
constant density solution); most of the others “describe 
situations which are frankly unphysical, and these do 
have a tendency to distract attention from the more use¬ 
ful ones.” [5]. One, the so-called Tolman VII solution 
appeared to have some physical relevance, but this was 
one of the solutions for which no explicit expression for 
the pressure was given. 


The Tolman VII solution has been rediscovered a num¬ 
ber of times and has appeared under different names, 
the Durgapal and the |Me hra| solutions being two 

examples. That these solutions can be used to describe 
realistic physical systems has been noted by many au¬ 
thors, including those of the two review papers mentioned 
above HO- It has been used as an exact analytic model 
for spherically symmetric stellar systems, and additional 
research has investigated its stability properties mm- 
While these later works were able to obtain the compli¬ 
cated expressions for the fluid pressure as a function of 
the radial coordinate, according to Finch and Skea [3] 
it still was not one of the “interesting solutions” since 
it lacked an explicit expression for the equation of state. 
The choice of parameters that has been taken by differ- 
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ent authors in order to completely specify the solution 
in many ways prevented the immediate interpretation of 
the physical conditions described by the solution. 

The reasons mentioned above are not sufficient to use 
or classify the Tolman VII solution as a physically viable 
one. Instead, we seek physical motivations for the via¬ 
bility of this solution, and indeed we find these in many 
forms: 


(i) From a Newtonian point of view, simple thermody¬ 
namic arguments yield polytropes of the form p(p) = 
kp 7 (here 7 is the adiabatic index sometimes written 
in terms of the polytropic index n, 7 = 1 + 1 /n, and 
k is known at the adiabatic constant that can vary 
from star to star) as viable models for neutron matter. 
When coupled with Newtonian hydrodynamic stability 
and gravitation, the result is the Lane-Emden differen¬ 
tial equation for the density profile, p(r). Solutions of 
the latter, obtained numerically, or in particular cases 
(7 = 00 , 2 , or 1 . 2 ) exactly, all have a distinctive density 
falloff from the center to the edge of the Newtonian 
star. This is a feature we wish physical solutions to 
have. Furthermore, this distinctive falloff is quadratic 
in the rescaled radius m, suggesting that even in the 
relativistic case, such a falloff would be a good first 
approximation to model realistic stars, which have a 
proper thermodynamic grounding. 


(ii) Looking at viable exact relativistic solutions to the 
Einstein equations, the one used extensively before 
1939 and even much later, was the Scliwarzschild in¬ 
terior solution. This solution has the feature that the 
density is constant throughout the sphere, and is not 
physical: the speed of sound (pressure) waves in its in¬ 
terior is infinite. However, this solution provides clear 
predictions about the maximum possible mass of rel¬ 
ativistic stars in the form of the Buchdahl limit m 
M < 4i?/9. The next best guess in this line of rea¬ 
soning of finding limiting values from exact solutions 
would be to find an exact solution with a density pro¬ 
file that decreases with increasing radius, since a sta¬ 
bility heuristic for stars demands that dp/dr < 0 , as 
expected from (i) in the Newtonian case. Extension to 
the relativistic Lane-Emden equation also requires m 
that (dp/dr)| r= o = 0, a property Tolman VII has. 


(iii) Additionally, an extensive review [T3j of most EOSs 
used from nuclear physics to model neutron stars con¬ 
cluded that a quadratic falloff in the density is a very 
close approximation to most such nuclear models-the 
differences of drastically different nuclear models from 
Tolman VII being only minor if only the density pro¬ 
files were compared. Since Tolman VII is precisely the 
unique exact solution to the full Einstein field equa¬ 
tions that exhibits a quadratic falloff in the density 
profile, we believe that it captures much of what nu¬ 
clear models have to say about the overall structure of 
relativistic stars. 


These three reasons taken together make a strong case 


for considering the Tolman VII solution as the best pos¬ 
sible exact solution that is capable of describing a wide 
class of EOSs for neutron stars. At the very least, it is 
as good a candidate that captures first-order effects in 
density of most nuclear model EOSs, and at best it is 
the model that all realistic nuclear models tend to, while 
including features like self-boundedness naturally, as we 
shall show. 

The purpose of this paper is to reexamine the Tol¬ 
man VII solution by introducing a set of constant param¬ 
eters that we believe provide a more intuitive understand¬ 
ing of the physical content of the solution. In addition, 
the solution now becomes a member of the set of “inter¬ 
esting solutions” since we provide an explicit expression 
for a class of equations of state derived from the solu¬ 
tion without any further assumptions about the matter, 
except for the Newtonian-like, and physically motivated 
quadratic falloff of the density. The EOSs will allow for 
further exploration of the predictions of the solution as 
well as a description of the material that makes up the 
star. The imposition of both the causality conditions 
where the speed of sound in the fluid never exceeds the 
speed of light and the different boundary conditions will 
provide further restrictions on the parameters associated 
with the solution. What this all leads to is a complete 
analytic model for compact stars that can be used to com¬ 
pare with recent observations of neutron star masses and 
radii. That the Tolman VII solution is consistent with all 
observations of astrophysical neutron stars leads to the 
conclusion that this exact solution is physically relevant 
while having features present in compact objects found 
in nature. 

This article is divided as follows: following a brief his¬ 
torical introduction in Sec.[TJ we re-derive the Tolman VII 
solution in Sec.[H] paying particular attention to the pres¬ 
sure expression with physically more intuitive variables. 
We then invert the density equation and use the pres¬ 
sure expression just found to derive a class of EOSs in 
Sec. Ill where we also carry out an analysis of the said 
class of EOSs. In the same section, we contrast the two 
different types of physical models that the solution ad¬ 
mits, and we shall also show how qualitative differences 
arise in the stars’ structure and quantitative ones appear 
in the predicted values of the adiabatic indices of the 
fluid. Finally we provide some brief concluding remarks 
in Sec. UV] 


II. THE TOLMAN VII SOLUTION 

Beginning with a line element in terms of standard 
areal (Schwarzschild) coordinates for a static and spher¬ 
ically symmetric metric, 

ds 2 = e^df 2 - e A(r) dr 2 - r 2 d(9 2 - r 2 sin 2 O&tp 2 7 ( 1 ) 

the Einstein equations for a perfect fluid source lead to 
three ordinary differential equations for the two metric 
variables v , A, and the two matter variables p and p. 
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However, these variables will not be the most practical 
ones to carry out our analysis. Instead two related met¬ 
ric functions, Z(r) = e -A( - r ) and Y(r) = e"M/ 2 , are in¬ 
troduced, as prescribed in Ivanov m- The reason for 


introducing these new metric variables is that with the 
assumption made for the density function, these variables 
will transform the original nonlinear differential equa¬ 
tions into linear ones which may then be easily solved. 


Given the metric equation ([I]) , the Einstein equations reduce to the following set of three coupled ordinary differential 
equations (ODEs) for the four variables Z,Y,p , and p: 


Kp = e 




+ 


1 

r 2 


Kp = e 




Kp = e 
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(2a) 

(2b) 

(2c) 


where the primes (') denote differentiation with respect 
to r, and k is equal to 87 t, since in what follows natural 
units where G = c = 1 are introduced. 


The first two equations (2a I and (2b) can be added 


together to generate the simpler equation 
, . 2ZdY Id Z 


( 3 ) 


which will be useful later on. In order to solve this set 
of ODEs, one begins with equation (2a) and assumes a 


specific functional form for the density, one that is mo¬ 


tivated from physical considerations according to (iii) 


Since this is a linear inhomogeneous ODE for Z(r), one 
can for the appropriately chosen form of p(r ) easily inte¬ 
grate this equation. The Tolman VII density has a simple 
functional form: 


P = Pc 


1 — P — 
n 


( 4 ) 


where the constant r b represents the boundary radius as 
mentioned previously, p c represents the central density 
at r = 0, and p is a “self-boundedness” dimensionless 
parameter that spans values between 0 and 1, so that 
when it is equal to 0, we have a sphere of constant density, 
and when it is equal to 1, we have a “natural” star, with 
density vanishing at the boundary. 

Although very simple, this quadratic function is known 
to provide a good approximation for the density pro¬ 
file of a number of neutron star’s EOSs. For example, 
Fig 5 in Ref. [13] plots the density profile of 12 EOSs and 
compares them to a function of the form given in equa¬ 
tion Q (for the p = 1 case). Therefore, the claim is that 
this functional form is a generic feature of many different 
types of nuclear EOSs and this suggests that at the very 
least some global features of such a density profile might 
describe the bulk properties of many compact objects. 

The set of three parameters that describe the density 
function will occur frequently in what follows and will 
be denoted as: II := {p c ,r b ,p}. The form of the density 


r 


function for p > 0 is physically realistic, since it is mono- 
tonically decreasing from the c ente r to the edge of the 
sphere, as argued previously in |(ii)| and (iii) in contrast 


to the constant-density exact solution frequently used to 
model such objects. 

Additionally, boundary conditions are required for the 
system, since we eventually want to match this interior 
solution to an external metric. Since the vacuum region 
is spherically symmetric and static, the only candidate by 
Birkhoff’s theorem is the Schwarzschild exterior solution. 
The Israel-Darmois junction conditions for this system 
can then be shown to be equivalent to the following two 
conditions m 


p(r b ) = 0, and 

Z(r b ) = 1 - — = Y 2 (r b ), 
r b 


(5a) 

(5b) 


where M = rn(r b ) is the total mass of the sphere as seen 
by an outside observer, and m(r) is the mass function 
defined by 


m(r) =47 t p(r)r 2 dr. 

J o 


( 6 ) 


Furthermore, the requirement of regularity for the mass 
function, that it must vanish at the origin of the ra¬ 
dial coordinate from physical considerations, leads to 
m(r = 0) = 0. On imposing (5b), one immediately writes 


Z in terms of the parameters appearing in the density as¬ 
sumption: 


y, x , /«Pc\ 2 , ( KPPc\ 

z(r) = 1 -(x) r + (w) 


=: 1 — br 2 + ar 4 . (7) 


In contrast, Tolman’s method was to assume the second 
form for Z (or equivalently for e _A ) in ([7]), and then 
obtain the density function from (2a) directly by differ¬ 


entiation. The physical constants p, p c , and r b occur fre¬ 
quently enough in the combinations shown above that 
the constants a and b as defined in Q will be used when 
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convenient. The solution methods for solving the ODEs 
obtained from the Einstein equations, particularly those 
leading to the Tolman VII solution, have been given in 
multiple references mm and will not be reproduced here. 

The complete Tolman VII solution is specified with the 
two functions ^ and Q below, together with the previ¬ 
ously given density function Q, and the metric function 
Z in equation ([T]): 

Y{£) =cicos(0f) + C2sin(0£), (8) 

where <j> = \J a/4. The quantity £ is a new radial variable 
whose explicit expression in terms of r is 


and it has been employed to simplify the expression of Y. 


Now that the full solution for the metric functions is 
known, the pressure can be computed through the re¬ 
lation (10), obtained from a simple rearrangement and 
variable change of ([3|: 


Kp(r) = 4 


VZdY 
~Y~ d£ 


Id Z 
r dr 


up. 


( 10 ) 


€(r) 



1 + \/l — br 2 + ar 4 


(9) 


This substitution results in a very complicated-looking 
expression for the pressure, 


np(r ) 


4</>[c 2 cos ((/>£) — ci sin (</>£)] \/l — br 2 + ar A 
Ci cos (</>£) + c 2 sin ((/>£) 


— 4ar 2 + 2 b — np c 



( 11 ) 


So far the two integration constants ci and c 2 appearing 
in the expression for Y, and therefore p, are completely 
arbitrary. Application of the the boundary conditions 
using equations (10), (5a), and (5b) leads to 




4j /Z&bf dV 

d£ 


-2 


£=£(> 


dZ 

da; 


, ( 12 ) 


where x := r 2 is another radial coordinate, and all the b- 
subscripted variables are the values at the boundary r = 
rt>. The cancellation shown results from matching to the 
exterior Schwarzschild solution. However, according to 
the second boundary condition (5a), the pressure has to 


vanish at the boundary; therefore equation (12) simplifies 
to 


i A Y 

K P l*=* t = 4 






which can be further simplified and rearranged as 


d Y 

b — axb 

np c 

(i 

A 

d£ 

«=«6 4 

4 

u 



(13) 


Since the ODE, equation (2c) for Y, is second order, a 


second condition is required. This is simply going to be 
condition (5b) restated as 


Y(x = x b ) = ^ 


1 — bxb + ax ? 


= l/l - K Pc r b 


=: 7 (14) 


The two equations (13) and (14) constitute the complete 


Cauchy’s boundary condition on Y. The integration con- 


stants ci and c 2 can now be determined from the simul- 

taneous equations: 


d Y 
d$ 

= cf> [c 2 cos (0£ b ) - ci sin (0&)] = a, 

£=£t 



.'. c 2 cos (</>£ b ) - ci sin (<j>£ b ) = a/<j>, 

(15) 

Y{£ = 

= Zb) = 7 



.-. c 2 sin ((/)£ b ) + ci cos (4>t;b) = 7- 

(16) 

This system can be solved by first multiplying (|15|) by 
cos (4>£,b), and (16) by sin(0£b), and adding the equa- 
tions obtained, yielding c 2 . Similarly, switching the mul¬ 
tiplicands and performing a subtraction instead yields Ci, 
and these can be given in the form, 


Q 

Cl = 7 cos (<fi£ b ) - - sin (^ b ), 

9 

(17) 


cx. 

c 2 = 7 sin (<j>£ b ) + - cos (0f&). 

9 

(18) 


The integration constants are ultimately computed in 
terms of the parameter set II, and in doing so this com¬ 
pletes the specification of the full Tolman VII solution in 
the new constant scheme. 

An important quantity to consider (since it establishes 
whether or not the solution is relativistically causal) is 
the adiabatic speed of sound waves in the fluid. The usual 
definition of this quantity in perfect fluids is v 2 = dp/dp. 
However, it will be convenient to find an expression for 
the sound speed directly from the differential equations, 
since the expression and functional form while completely 
equivalent is simpler to work with. First, from the ex¬ 
pression for the density Q, one obtains the derivative 

dp 2 pp c 

~r = —2 a 

dr r / 
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which is zero only at r = 0. For the other equation, the 
conservation of the energy-momentum tensor V {T l j = 0 
reduces to 


dp _ v'(jp + p) _ (p + p) dY 
dr 2 Y dr ’ 


(19) 


in the j = 0 case. These two expressions can be used to 
find dp/dp for every value of r but the center, so that 


y 2 = dp = 

dp V dr / dr 


dp /dp\ r 2 (p + p)dY 
2 pp c rY dr 


( 20 ) 


Using the expressions for all the terms in this formula, 
one obtains a relatively simple expression for the speed 
of sound. 

The bulk modulus K of a fluid is a measure of the 
resistance of a fluid to change its volume under an applied 
pressure. For perfect fluids it is related to the speed 
of sound in the media through I\ = pv 2 . This is also 
a quantity which may be computed for the fluid in the 
interior. This calculation demonstrates that the order 
of magnitude of the bulk modulus is significantly higher 
than any currently known substance by many orders of 
magnitude. 

The next step to understanding this solution is to in¬ 
vestigate the behavior of the solution as the parameter 
set is varied. The particular choice of parameters will 
be those that are associated with what one might expect 
for realistic compact astrophysical objects. As a result, 
central densities p c ~ 10 15 g • cm -3 will be typical. Simi¬ 
larly, radii rt ~ 10 6 cm (i.e. 10 km) will often be used for 
the same reason. As stated above, the density profile Q 
will decrease quadratically, and this provides a good ap¬ 
proximation of what one would expect from a number of 
neutron star EOSs. Figure [T] plots the density as a func¬ 
tion of radius for different values of p. which controls the 
relation of the surface density to the central density. 

The surface density ranges from a zero value when 
p = 1 to increasingly higher densities as p is decreased. 
In the literature m, models having zero surface densities 
have been called “natural,” and those with nonvanishing 
surface densities have been called “self-bound.” It is for 
this reason that p is called the “self-boundedness” pa¬ 
rameter. 

Similarly, the complicated expression for the pressure 
given by equations ( |TT] ), and p~8j ) can also be plot¬ 
ted as a function of the radius. Of importance here is 
the fact that while the densities might not vanish at the 
boundary rt,, the pressure for all parameter values must 
do so according to the boundary condition |5a| ). This 
is eminently clear in Fig. [2] where we see the pressures 
associated with the density curves shown in Fig. [T] Simi¬ 
larly the speed of sound and bulk modulus, all associated 
with the matter content in the star are shown in Figs. [3] 
and [4] respectively. 

The functions Z(r) and Y(r) representing the solu¬ 
tions to the differential equations (2) are given in Figs. [6] 
and [5] respectively, again for different values of the self¬ 
boundedness p. Equivalently the metric coefficients in 



Figure 1 . Variation of density with the radial coordi¬ 
nate inside the star. The parameter values are p c = 
1 x 10 15 g • cm -3 , rt = 1 x 10 6 cm and 0.6 < p < 1.0 . 



Figure 2 . Variation of pressure with the radial coordi¬ 
nate inside the star. The parameter values are p c = 
f x 10 15 g • cm -3 , rt = 1 x 10 6 cm and 0.6 < p < 1.0 . 



r(m) 


Figure 3 . Variation of speed of sound with the radial co¬ 
ordinate inside the star. The parameter values are p c = 
1 x 10 15 g • cm -3 , rt = 1 x 10 6 cm and 0.6 < p < 1.0 . 


Schwarzschild form: the form most often used in the lit¬ 
erature for specifying static spherically symmetric mod¬ 
els can be obtained from Y(r) and Z(r). For the sake of 
completeness, A(r) is plotted in Fig.[7]and v(r) is plotted 
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Figure 4 . Variation of bulk modulus with the radial co¬ 
ordinate inside the star. The parameter values are p c = 
1 x 10 15 g • cm -3 , rt = 1 x 10®cm and 0.6 < p < 1.0 . 



Figure 7 . Variation of A metric variable with the radial co¬ 
ordinate inside the star. The parameter values are p c = 
1 x 10 15 g • m -3 rt = 1 x 10 6 cm and 0.6 < p < 1.0 . 


in Fig. [8] 



Figure 5 . Variation of Y metric variable with the radial 
coordinate inside the star. The parameter values are p c = 
1 x 10 15 g • cm -3 , rb = 1 x 10 6 cm and 0.6 < p < 1.0 . taking 
the various values shown in the legend. 



r (m) 

Figure 6. Variation of Z metric variable with the radial co¬ 
ordinate inside the star. The parameter values are p c = 
1 x 10 15 g • cm -3 , Vb = lx 10 6 m and 0.6 < p < 1.0 . 


u 

OJ 

E 



r (m) 


Figure 8. Variation of v metric variable with the radial co¬ 
ordinate inside the star. The parameter values are p c = 
1 x 10 15 g • cm -3 rb = lx 10 6 m and 0.6 < p < 1.0 . 


The redshift z s of light emanating from a star as per¬ 
ceived by distant observers is another quantity that po¬ 
tentially can be measured. This quantity can also be 
calculated in our model, from the relation 





2 m{r b ) 
n 


2 


- 1 . 


The redshift value at the surface of the star for different 
values of /.t is shown in Fig. [9j 


III. THE EQUATION OF STATE AND 
PHYSICAL MODELS 

A nice feature of the density assumption 0 is that it 
can be inverted to easily obtain r as a function of p. This 
allows one to generate an equation of state (EOS) for this 
solution. The full equation of state is given below: 

P(p) = ~2 oJ tl fo 2 | /l1 “ h ' 2 \!■ft c°t 2 f 2 + 47 r/ii/i 2 p| , 
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Figure 9 . The redshift z s at the surface of the sphere for 
different values of p. 


where fi(p) and f 2 (p) are functions of the density: 

/i(p) = 50 — 3 - ~J~P + 327r 2 /i?p 2 

and 

VS/i/12 + /ii — 167r/ii/i 2 p 
20 /i 2 C 

The constants hi and /i 2 are determined by the central 
density and p, as follows: 


f2(p) = 


hi 


n 


2-n-pcP 


and 


h 2 = 


3 

8 t rpc 


while the constant C is expressible as a complicated func¬ 
tion of the parameters only, in terms of the auxiliary 
variables a and \\ 


C = 



I /ii(4fe 2 -/ii) 

8 rlh 2 ~hf + x 


exp 


arctan 



with 


X = 4 \]h 2 {Ah 2 r^ - h\rl + h\h 2 ), 
a = 16 h 2 r\ + 8Trp c h 2 h 2 (l — p) — 2 h\. 


It should be noted here that no assumption about the 
nature of matter, except for the very general thermody¬ 
namic prescription of a perfect fluid, has gone into this 
solution. Everything else, and in particular the equation 
of state, was obtained solely by virtue of the held equa¬ 
tions and the density profile Q. With the equation of 
state now given explicitly as p = p{p), it is a simple mat¬ 
ter to find the derivative dp/dp for the speed of pressure 
waves, and this yields precisely the same function as the 
one found previously in equation (20). 

The expression for this class of EOS is somewhat com¬ 
plicated, but it is not without physical interpretation, 
contrary to what Tolman [3] thought in 1939: 


The dependence of p on r , with e~ x / 2 and e~ 1 ' 
explicitly expressed in terms of r, is so com¬ 
plicated that the solution is not a convenient 
one for physical considerations. 


By virtue of having a class of exact EOS, there is the 
possibility of two separate interpretations for an EOS 
that arise from the analytic expressions. This classi¬ 
fication can be seen as a practical way of interpret¬ 
ing a class of EOS that has four different parameters, 
not all independent of each other. Both p(p: II) for 
Pb = Pc{ 1 — At) < p < Pci with the values of the elements of 
n, in particular p c , fixed (henceforth called EOS1); and 
p{p = p c \ II), with the parameters of II varying between 
limits imposed by causality (EOS2), could be candidates 
of the EOS. 

In the literature, both interpretations have been used, 
and sometimes even interchanged. However, each has a 
completely different content in that the first interpreta¬ 
tion expresses how the pressure of the fluid changes in 
moving from the center of the star r = 0, p = p c , to the 
boundary r = rb,Pb = Pc{ 1 — p), while all the integra¬ 
tion constants, and hence parameters II, are kept fixed. 
Those seeking an interpretation of a unique EOS that 
should be applicable to all neutron stars without excep¬ 
tion would find this interpretation sufficient. 

The second interpretation, by contrast, looks closely at 
the fluid material itself and how the pressure at a certain 
point in the star changes as the density of the fluid at 
the center changes. Given that the central mass density 
of a compact star is inaccessible, this interpretation is 
of interest to those who believe that the central density 
should be a free parameter in a neutron star model. This 
would allow one to explore the possibilities that such a 
parameter change has on the observable quantities of the 
stars. 

At this point in the derivation, a causality condition 
has not been imposed upon expressions of EOS1 or EOS2. 
Therefore, the class of EOS obtained can take a wide 
range of parameter values, as long as the metric functions 
and derived curvature tensors do not have singularities. 
The conditions leading to such singularities are narrower 
than the causality criterion, and enforcing the latter en¬ 
sures that the parameter values do not cause singular 
behavior in the solution, and hence the class of EOS. 

We first carry out an analysis of EOS1, and find that, 
to a high degree of accuracy, the variation of p{p\ II), with 
p, and equivalently r, is very close to that of a generalized 
polytrope of the general form p = kp 7 — po , where po is 
a pressure constant chosen such that p vanishes when 
p = pb at the boundary of the star and, as usual,7 and 
k are the adiabatic index and the adiabatic constant, 
respectively. This relation is very obvious from the shape 
of the curve in the “natural” p = 1 case as is seen in 
the one curve in Fig. [TOJ and all the curves in Figs. [T2] 
and 13 It is interesting to note that the p < 1 cases 
all show a behavior similar to that found for other self¬ 


bound EOSs. (See e.g. Fig. 1 in the review by Lattimer 
and Prakash| [IB] ). Indeed, Fig. 13 even seems to suggest 






































that varying the boundary radius rb changes the value of 
k in the polytrope, and Fig. 12 that varying p c changes 
the value of 7 in the polytrope. 



Energy Density (g/cm 3 ) 


Figure 10 . Log-log plot of pressure versus density for neutron 
star models determined by different p, but the same p c and 
tv The densities and pressures are in cgs units, and the II is 
fixed by the following: rb = 10 6 cm, p c = 10 15 g • cm -3 . Since 
pressure is a decreasing function of distance from the center, 
large densities indicate points closer to the center of the star. 

Models employing polytropic perfect fluids use similar 
values for the adiabatic index 7 as what we find for a 
range of different values of parameters II. This is shown 
in Fig. ITT] which treats 7 as a continuous variable defined 
by 7 = dp°gp] , and can be understood as the slope of the 
previous log-log graph. 



Figure 11 . The adiabatic index variation from the center to 
the boundary of the star for different values of the parameter 
p. The other parameter values are the same as those in Fig. 110 | 

From this figure it becomes evident how both types 
of stars have an interior structure well described by a 
polytrope with an adiabatic index close to 2.5. The “self¬ 
bound” stars exhibit the existence of an envelope consist¬ 
ing of material that is considerably stiffer than that found 
in the interior. Physically this is intuitive: for fixed p c 
and 77, the self-bound stars will become more and more 
massive as p decreases. The increasing boundary density 
discontinuity requires a stiffer exterior mass distribution 


to maintain the hydrostatic equilibrium condition. 

A notable characteristic of the class of EOS is its un¬ 
canny ability to distinguish between the different types 
of matter that make up the natural and the self-bound 
stars. Since the mass density Q is a monotonically de¬ 
creasing function of stellar radius, Fig. [Tl] can be thought 
of as the equivalent of a “flipped” and “rescaled” plot of 
the adiabatic index as a function of density. For the case 
p = 1, Figs. 10 and 11 are consistent with a number of 
hadronic EOS proposals. For densities in the range of 
1-10 xl0 14 g-cm 3 the adiabatic index ranges from 2.7 
for the most dense nuclear material to 2.0 for the lower- 
density material. This type of behavior is found, for 
example, in a model proposed byGlendenningjT6j, that 
consists of a mixture of baryons. 

For the self-bound models, where p < 1, the Tolman 
EOS is consistent with quark models using the MIT bag 
model. Fig. [To] is similar in nature to the strange quark 
models (SQM1-3) shown in Fig. 1 in [T3] while Fig. IT 
has similarities to the work of |Casali and Menezes[ who 


analyze the MIT bag model (see also the book by Haensel 


et al. PU) and where it is found that the EOS for quark 


matter stiffens significantly at low densities jT7]. The 
adiabatic index of that material reaches very high values 
(e.g. 7 > 7 for densities less than 10 14 g • cm -3 ). 

To fully understand the nature of the EOS, the ef¬ 
fect on the matter resulting from changing the other two 
parameters of II can be investigated. Figure 12 demon¬ 
strates how the EOS inside the star changes as the central 
density p c changes, and similarly Fig. [I3| shows how the 
EOS varies with changes in the magnitude of the bound¬ 
ary radius r&. 



Figure 12 . The effect of changing the central density p c 
through 1 order of magnitude on the EOS with a boundary 
radius kept fixed at 10 km, for the natural EOS with p = 1. 
We notice that the slopes of the lines change by very little. 
The legend also provides the corresponding masses associated 
with each parameter choice in solar mass units. 

What these figures show is that the value of the free 
parameters can change by up to an order of magnitude, 
yielding drastically different masses, while maintaining 
the same general polytropic behavior, independently of 
the parameter choices. At the highest densities, Fig. [12] 
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Figure 13 . The effect of changing the boundary radius Tb 
through 1 order of magnitude, while the central density is 
kept fixed at 1.5 x 10 15 g-cm -3 . We notice that the lines 
remain mostly parallel, suggesting that only the k value in the 
polytrope is changing, while the adiabatic index is remaining 
about the same. The masses of the corresponding stars are 
also provided in solar units in the legend. 


shows that the EOS are nearly independent of the central 
density parameter. Figure [13] which demonstrates how 
the EOS changes under changes in the boundary radius, 
indicates that the adiabatic index 7 is fixed over a large 
range of densities, but as expected the adiabatic constant 
k is different for different stars. The self-boundedness pa¬ 
rameter p, however, changes the character of this poly¬ 
trope very much, as is clear in Fig. ED hinting that the 
same EOS can have a richer structure than can be ulti¬ 
mately specified by the central density alone. 

This polytropic behavior is very satisfying, since we 
started by trying to model a relativistic star from a New¬ 
tonian picture. That a class of EOS globally similar to 
the solutions of the Lane-Emden equations becomes ap¬ 
parent when we extend the nonphysical Schwarzschild 
interior to a more realistic density profile suggests that 
Tolman VII is at least as good as the Newtonian neutron 
stars, however with relativity being taken into account. 

Now turning to the second way to characterize the class 
of EOS, concentrating on the behavior of the fluid mate¬ 
rial itself, independent of the geometry of the star, we de¬ 
termine how different physical quantities depend on the 
values of the central density p c . The total mass-energy is 
defined as 


M = 4-7T 



4tt r\p c (5 - 3/i) 
15 


( 21 ) 


The mass is important, since it is the only directly and 
reliably measurable quantity that can be obtained from 
neutron star observations. Lattimer and Prakaslr|l31 fl9 l 
[2D] and others EUI221 have ruled out certain EOS2 based 
on mass and spin measurement of neutron stars. The for¬ 
mer have also used Tolman VII, to constrain other EOS2 
based on nuclear microphysics, and have even postulated 
that Tolman VII could be used as a guideline discrimi¬ 
nating between viable and nonviable EOS2 [20]. If this 


postulate is true, given that the complete Tolman VII 
EOS2 is known, the condition that the solution must be 
causal can be applied, independent of measurements first, 
and then compared with the previous works [20] [22]. 

This is done in Fig. |14[ where we superimpose the re¬ 
sult of Ref. [22], on our own analysis of the whole solution 
space II. The surface shown is that for the values at which 
the speed of sound v s = {^/Ap/Ap 

the fluid sphere just reaches the speed of light. This is a 
sufficient condition for the solution to be causal, since v s 
is a monotonically decreasing function of r in the sphere. 
Any point located below this surface has coordinate val¬ 
ues for M, p Cl and p that represent a valid causal solution 
to the Tolman VII differential equations. The orange line 
is the previous result obtained by Glendenning l22j from 
rotational considerations. 


at the center of 

r= 


Causal surface 


3 M© 



Figure 14 . The mass of possible stars just obeying causality. 
The grey surface obeys the equation v s (r = 0 ) = c. Every 
point below the surface is a possible realization of a star, and 
we can potentially read off the mass, central density, and p 
value of that star. The numbered lines represent stars with 
the same mass that are causal, i.e. they are projections of the 
causal surface onto the p c -p plane. Glendenning’s [ 22 ] curve 
is shown in orange and represents a limit in the natural case 
only, and according to our results is acausal, being above our 
surface. The p = 1 plane’s intersection with our graph is the 
graph given in Ref. i 2 U], and here too our prediction is more 
restrictive. 

Imposing causality to constrain the parameter space 
n is not a new idea. However, having an explicit EOS 
allows one to easily generate the causal surface shown 
in Fig. |14[ exactly without having to do any numerical 
gymnastics to find the speed of sound. 

One method often used to distinguish between differ¬ 
ent EOS2 has been to calculate the compactness ratio, 
given by (3 = . In the case of Tolman VII the values 

of /? for a large range of parameter variations II are rela¬ 
tively constant. This means that even though we might 
change the value for II of the stars, the ones bordering on 
causality share very similar compactness, albeit one that 
is lower than that previously thought possible. We show 
how this compactness (3 varies with p in Fig. |15[ Clearly 
the variation of (3 is small (0.27 < (3 < 0.29) over a large 
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Figure 15 . The compactness as a function of the self¬ 
boundedness parameter /x. This plot was generated by vary¬ 
ing rb from 4 km to 20 km for fixed /x and finding p c and 
subsequently the compactness each time, such that the sound 
speed was causal at the center of the star. The curve shown 
is a polynomial fit, and the box-and-whisker plots(very small 
in green) show the variation of /? for fixed /x, but different rb. 
The very small whiskers justify the pertinence of /3 as a useful 
measure in the analysis of the behavior of the model. 



r b (km) 


Figure 16 . The mass M in solar units versus radius r b in 
kilometers of a few stars for which these values have been 
measured. We use error bars to denote observational uncer¬ 
tainties, and colored bands in the case where only the mass 
is known. The lines we show and the ones that are at the 
causality limit in the Tolman VII model, and the whole area 
below those lines can be causal parameter choices for stars. 


range of values for p (0.45 < p < 1.0). A value for the 
maximal compactness of about 0.34 from rotational and 
causality criteria was obtained by Ref. [20]. Our analysis 
shows that /3 is less than 0.3 for all possible stars, if Tol¬ 
man VII is a valid physical model (at least as a limit) for 
stars. Recently radius measurements of a limited num¬ 
ber of neutron stars have been obtained [281426] . These 
are shown along with some other stars of known mass in 
Fig. [Tg] We also superimpose a few of the limiting causal 


curves obtained for different values of p from Tolman VII, 
to show that Tolman VII is not ruled out by observational 
results, even though it predicts lower compactness than 
most nuclear models. The dotted lines shown in the fig¬ 
ure represent the causal limits for different values of the 
self-boundedness parameter /x. 

For a fixed value of /i, causality requirements deter¬ 
mine the relationship between the central density and the 
boundary of the star. A typical M-R curve for a particu¬ 
lar EOS2 would be determined by fixing the parameters 
associated with a particular EOS2 and then computing 
the mass as a function of the star’s radius. All possible 
curves of this nature lie to the right of the curves shown 
in Fig. [TG] and thus any curves in that region represent 
a viable model for neutron star structure obtained from 
a Tolman VII solution. This is one advantage of hav¬ 
ing a complete analytic solution to the Einstein equa¬ 
tions. Rather than integrate (using numerical methods) 
the TOV equation for a specific set of EOS parameters, 
an entire class of solutions are provided by the complete 
analytic expressions. Physical constraints on the solu¬ 
tion (such as causality) then provide restrictions on that 
class of solutions. What the curves shown in this plot 
represent are the maximum possible cases for the (M, r^) 
relation provided by the causal Tolman VII EOS solu¬ 
tions. Therefore, the whole space to the right of those 
lines can potentially yield a causal EOS and thus Fig. [TG] 
shows that Tolman VII is not ruled out yet as a viable 
model for compact objects unless the radii of some of the 
neutron stars with known masses are so small that they 
would be represented by a point to the left of the curves. 

The lines shown are on the edge of causality in the 
following way: the speed of sound, a monotonically de¬ 
creasing function of the radius, is just equal to the speed 
of light at the center where r = 0- that is, the lines are 
the counterparts of those points that make up the sur¬ 
face of Fig. [14] Since all observations of compactness 
are bounded by the most extreme Tolman VII model, we 
claim that the solution may actually realized by compact 
stars in nature. 


IV. CONCLUSION 

A complete analysis of the Tolman VII solution has 
been carried out, and it was found that it is a physically 
valid solution with a huge potential for modeling physi¬ 
cal objects. The class of EOS this solution predicts has 
been found, and in certain regimes behaves very much 
like a polytrope with an adiabatic index of 2.5, inde¬ 
pendently of the choices of two parameters: the central 
density and the radius of the star. The third parameter, 
the self-boundedness, changes the polytropic index dras¬ 
tically, particularly at the edge of the star, as expected 
from a naive Newtonian approach to stellar structure. 
That this solution has a density profile that is very close 
to the hand-picked, but thermodynamically motivated 
polytropes of Newtonian stars, while still being a full 
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relativistic model is a very good reason to take its pre¬ 
dictions of stellar structure seriously. 

It is also interesting that the type of matter that pro¬ 
duces the different EOSs depends crucially on the value of 
p. Hadronic matter is obtained with fi = 1, while /i < 1 
stars would appear to be made up of quark matter. 

Using the EOS derived from Tolman VII, we are able 
to compute the speed of pressure waves, and imposing 
causality on the latter results in a more restrictive limit 
on the maximum compactness of fluid spheres allowable 
by classical general relativity. This is possible to do with¬ 
out the use of numerical computations because of the ex¬ 
act form of the class of EOSs generated by the quadratic 


density falloff assumption in Tolman VII. 

The solution is, moreover, stable under radial pertur¬ 
bations, since the speed of these pressure waves is finite 
and monotonically decreasing from the center outwards, 
thus satisfying the stability criterion in Ref. [27] . If we 
believe as in Ref. [20] that Tolman VII is an upper limit 
on the possible energy density p c , for a given mass M. 
some known models m which predict higher compact¬ 
ness than Tolman VII will have to be reconsidered, since 
they still maintain a quadratic density profile to first or¬ 
der, and thus cannot also be causal at these higher com¬ 
pactnesses. 
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